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^ A more  genera]  mathematical  model  has  been  developed  for  the  prediction  of  signal  power 
generated  by  the  Coherent  Anti-Stokes  Raman  Scattering  process  (CARS).  Like  previous  treatments 
the  model  assumes  that  the  laser  beams  have  Gaussian  profiles.  However,  provision  is  made  for 
varying  a number  of  experimental  parameters,  either  singly  or  in  combinations.  In  particular,  the 
model  allows  for  non-collinearity  of  the  laser  beams,  different  beam  diameters,  lateral  offset  of  the 
beam  centers,  and  relative  displacement  of  the  focal  points  of  the  beams.  .. 
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20.  Abstract  (Continued) 

Computational  time  requirement!  of  the  new  model  an  not  significantly  greater  than  pmvious 
treatments  of  identical,  colUnear  Gaussian  beams,  although  the  numerical  integrations  must  be 
handled  more  carefully.  Significant  predictions  of  the  model  are: 

1.  Substantial  power  gainscaft  be  realized  by  tightening  the  focus  of  the  (dye  laser)  beam 
relative  to  the  primary  pumping  beam. 


2.  Signal  intensity  is  enhanced  by  a slight  displacement  of  the  focus  of  the  w2  beam  in  the  direction 
of  propagation. 

3.  Substantial  deviations  from  colllnearlty  decrease  the  signal  intensity  sharply,  but  a small  offset 
results  in  a slight  enhancement. 


4.  Crossing  the  laser  beams  decreases  the  signal  power,  but  the  magbitude  of  the  effect  is  smaller 
than  might  be  expected. 

5.  The  loss  of  signal  power  which  results  from  a non-zero  crossing  angle  can  be  compensated  for 
(up  to  a crossing  angle  of  about  1°)  by  introducing  a small  lateral  offset. 

6.  For  a typical  set  of  experimental  parameters,  approximately  80%  of  the  CARS  power  is  generated 
in  a region  1 cm  long  with  its  center  at  the  beam  focus. 
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MODEL  FOR  PREDICTING  RELATIVE  SIGNAL  POWER  IN 
COHERENT  ANTI-STOKES  RAMAN  SPECTROSCOPY 


I . INTRODUCTION 

With  the  current  availability  of  versatile  laser  sys- 
tems, there  is  increasing  interest  in  higher-order  Raman 
scattering  techniques  as  tools  for  studying  a wide  variety 
of  chemical  systems  (1) . These  methods  have  numerous  ad- 
vantages over  conventional  Raman  spectroscopy,  such  as  in- 
creased spatial  resolution,  greatly  enhanced  signal  inten- 
sities, spatial  filtering,  etc.  (2). 

Typical  of  the  higher-order  Raman  processes  and  one  of 
the  most  promising  from  the  standpoint  of  practical  applica- 
tions is  the  so-called  coherent  anti-Stokes  Raman  spectro- 
scopy (CARS) . In  this  case  a sample  is  irradiated  with  two 
laser  beams  of  frequency  wJ,  and  , and  a signal  corres- 
ponding to  a characteristic  Raman  transition  of  the  sample 
is  generated  at  the  frequency  ol^  = 2u/,  - Uj\  . The  signal 
is  coherent  and  many  times  more  intense  than  the  correspond- 
ing output  obtained  by  conventional  Raman  spectroscopy  (3) . 
An  additional  advantage  is  realized  when  the  technique  is 
applied  to  condensed  phases  in  which  there  is  a significant 
amount  of  dispersion.  The  phase-matching  conditions  of  the 

Note:  Manuscript  submitted  February  8, 1979. 


multi-photon  process  require  that  the  two  laser  beams  cross 
at  a certain  angle  in  the  active  region  of  the  sample.  The 
result  is  that  the  coherent  signal  emerges  in  a direction 
different  from  that  of  either  of  the  two  laser  beams,  and 
thus  can  be  separated  without  the  use  of  a monochromator. 
For  gases  the  optimum  crossing  angle  is  near  zero,  so  that 
spatial  filtering  can  be  achieved  only  by  sacrificing  sig- 
nal intensity. 

Theoretical  calculations  of  CARS  signal  strengths  in 
gaseous  media  have  yielded  very  optimistic  predictions  of 
the  ultimate  sensitivity  of  the  technique.  The  first  rig- 
orous treatment  of  the  conversion  process  for  collinear, 
focused  Gaussian  beams  was  carried  out  by  Bjorklund  (4) . 

His  calculations  yielded  results  in  general  agreement  with 
the  experimental  observations  under  the  conditions  of  the 
mathematical  model . 

For  realistic  predictions  of  the  optimum  CARS  perfor- 
mance under  actual  experimental  conditions,  a more  general 
mathematical  model  is  required.  In  a typical  experiment 
the  laser  beams  are  not  precisely  collinear  and  do  not  have 
identical  profiles.  Furthermore,  the  beams  may  be  deliber- 
ately crossed  in  a gaseous  sample  in  order  to  achieve  in- 
creased spatial  resolution  in  the  active  region  and  spatial 
separation  of  the  signal  at  the  detector.  It  is  important 
to  know  the  extent  to  which  performance  is  degraded  when 
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the  experimental  parameters  are  varied  away  from  the  ideal 
case.  Therefore,  a more  flexible  predictive  model  would  be 
of  substantial  value  in  selecting  the  optimum  experimental 
conditions . 

In  this  work  the  treatment  of  Bjorklund  has  been  ex- 
tended to  include  the  effects  of  non-collinear , non-confo- 
cal  beams  with  different  profiles.  In  principle,  the  model 
is  capable  of  dealing  with  arbitrarily  large  crossing  angles 
between  the  two  laser  beams,  but  it  has  been  applied  only 
to  the  range  of  practical  interest  in  gases  (Sc  £ 3° ) . 
Modifying  the  model  to  allow  for  the  treatment  of  larger 
crossing  angles  would  require  additional  algebraic  manipu- 
lations, but  would  not  increase  the  computational  time  sig- 
nificantly. 

II.  THE  THEORETICAL  MODEL 

A.  The  General  Case 

In  this  section  the  most  general  experimental  arrange- 
ment covered  by  the  model  is  defined.  All  of  the  special 
cases  considered  below  will  be  obtained  from  the  general 
model  by  appropriate  restrictions  of  the  parameters. 

Let  us  start  with  a laser  beam  of  frequency  Ul  ( with 
propagation  vector  directed  along  the  longitudinal  axis  of 
the  sample  cell  (taken  to  be  the  z-axis) . Assume  the  beam 
has  a Gaussian  profile  (5) , and  assume  it  is  focused  inside 
the  cell  so  that  the  electric  field  strength  is  given  by 
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E , < r ) - 


- ll-ti— ) exp  ~ j 

b‘;  rL  t,0+^./b.)J 


where 


b.  = i,r,i 


is  a reduced  form  of  the  radius  f,  of  the  beam  waist  and 

€>(*)  = l(Z-f i)  . (3) 

The  parameter  "f|  is  the  z-coordinate  of  the  beam  waist. 
Without  loss  of  generality,  the  origin  of  the  coordinate 
system  can  be  taken  to  be  at  the  focus  of  the  beam,  so  that 

Z{(l)=Zi  . (4) 

Let  the  second  laser  beam  have  frequency  id.,  and  propa- 
gation vector  . The  beam  crossing  geometry  is  shown  in 
Figure  1.  Note  that  the  x-axis  is  chosen  in  such  a way  as 
to  make  both  fl  f and  parallel  to  the  xz-plane.  The  y^ 
axis  is  chosen  to  be  parallel  to  the  y-axis.  The  xxzi~ 
plane  cuts  the  y-axis  at  y = yc  , and  the  focus  of  the  sec- 
ond beam  (zz  = 0)  is  displaced  from  the  "crossing  point" 

(y  = yc,  z = zt)  by  the  distance  d. 

Assume  the  second  beam  is  also  a focused  Gaussian.  In 
the  coordinate  system  (.xx  iYx  >zx)  its  electric  field  vector 


is  given  by 


where  OL  is  an  arbitrary  phase  factor  and  all  other 
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quantities  have  meanings  which  are  analogous  to  the  corres- 
ponding parameters  in  Eq.  (1).  For  example, 

£l(l  = (6) 

where  is  the  distance  between  a particular  point  in 
space  and  the  focus  of  the  second  beam,  measured  in  the  di- 

.Jk 

rection  of  The  focus  of  the  second  beam  is  assumed  to 

lie  in  the  plane  zx  = 0. 

By  making  use  of  simple  geometry,  it  is  possible  to 
express  E2(r)  in  terms  of  the  (x,y,z)  coordinate  system. 

The  equations  of  the  coordinate  transformation  are 


xx  = xcos3c  - (z  - zc)sin5e 

(7a) 

0 
>1 

1 

>1 

II 

>1 

(7b) 

z2  = xsin3c  + (z  - zc)cos9c  - d 

, (7c) 

and  the  result  is 

(:*(?)  = £i0 

(l+iEi/l»zy 

Ax  C 

] . <8> 

where 

tzU)  - + Ki'ic)  co<>  -Id, 

(9) 

and  bx  - is  the  reduced  radius  of  the  iA)j.  beam. 

If  it  is  assumed  that  both  laser  beams  are  polarized 
along  the  y-axis,  then  the  polarization  of  the  medium  at 
the  frequency  1/0^  = Z^\-t)jXxs  given  by  (6) 


P.,  if)  = \ NX (Zuii-iOj.)  £*0  Ej,0  e”1*  bf  bz 

*C*p[”l^*s,’rt^c  + iC*sk)]{b|  +t,£i)'i(bx-‘  £*)'  ' 

> d^L-  z*«  ft1*  ua)Ab,-»i  6i)] 

* «-X.jo  j-J|— + ] > tl0) 

where  N is  the  number  density  and  y,(Iu)\  -u)*)  is  the 
third-order  polarizability  of  the  medium  at  the  frequency 


2.uj(-  The  Fourier  transform  of  Pj(r)  is 

f l&)-(  iff)  •[dl ' P,  ( x ')  £ 1 * ' r 

-Co  -CJ  -CD 


(ID 


Assume 


Kx,  « K2  and  |K|  - k5  « (K| , k,  , (12) 

so  that  Kleinman's  solution  (7)  to  the  wave  equation  is 
valid.  These  conditions  are  well-satisfied  for  laser  beams 
at  optical  frequencies.  Therefore,  the  Fourier  components 
of  the  radiation  generated  in  the  active  medium  by  the  CARS 
process  are  given  by 


7-TT  i l 


k)<1K  r i 


(13) 


where  fto  is  the  vacuum  propagation  number  of  the  generated 
radiation  and 

, (‘  -dip 

(Oi)Ac  r 


t 


(14) 


The  total  generated  field  at  the  frequency  u)j  is 


*~V> 
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obtained  by  substituting  the  expression  for  P CK) , Eq.  til) , 
into  Eq.  (13)  , using  Eq.  (14) , and  integrating  over  K: 

,ct  /CO  roo 


( r ) = ^,^x(^K1(dK2  E3-  (r ) 

-Cl  -00  --OP 

Defining 

A#  - \ N X C iuj » -u)z)  E*0  E1q  t lotb,  b2 

we  obtain 

,00  /.CD  r| 


(15) 


(16) 


-oo  -oo  o 

-*  * r<D  /'  ■*  / CO 


-Oo  -'.oo  -00 

'Ij'i 


X ^i>,  + i£,U')j  {b2- i ti  Cx\Z')\ 

-iK-r 

b b,  + ^ £ , £2')  J 


-1 


x ex' 


P 


u 


[x'CD-i^c  -(Z'-2c)5i«3c  ]* 


(17) 


The  first  step  is  to  integrate  over  Kf,  using  the  re- 
lationship (8) 
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to  reduce  the  integral  to  the  form 
; < * i roc  r rl 


-to  ■'-oa 

v/><r  /■■*>  r0^' 


^ht*  i dy\di‘  £ 

-CO  J-00  \oj 


yc  x* 


The  integration  over  z'  can  then  be  carried  out  to  yield 
E 5(f) - & A' 0 


-co  J-a? 


x«i>[; di,-iitn$cX i-f) [bi-itji'.d-fti 

r z i, cvjri  ] r **■  fr 

* °r  L-  rauT^ii  H r w-.tav,o-^H[x  t05Sc 


Finally,  the  integration  over  y'  is  carried  out, 


(22) 

also  using 


Eq.  (21) . The  rather  complicated  expression  which  results 
can  be  simplified  by  defining 


G,  Ip)  = b|  +i  £,Ui-p)iJ  “ + 2i(i-p)2  (23 


and 


(24 


Then 


E3(  r)'  (7)  Nfc*  t J 

> ^ -2^,+^zco!.9c)i|o 


S'  **  • A . 

, -L  E.i<i»3c  X*  . 

e Go.i/.p) 


eKp (- ' 


-Cc 


f 


p.  z-d.  U3v  ^jr  Hi  1 

' eY U ) + *2f  '*  JL"6,tr)"6t<*-,ft  “pj 


.1/ 


T 


i dx 


fi. filial  L^*  -iiT1" 

iG-jCv'.f)  + “ip" I l G,tf)  GiCx.f)  lip 


Let 


D liftp 

D|tp)  - - — t +-  77  J 

' G|(  j>)  ** 


(26) 


so  that 


E*<t)=(:)VA'eU2*,"Wlei^‘-a‘fc,i 

* (>"'  Gr.tp1  + ‘ +tt)ip 

^ [d)£,«Kp[-Li^X^5c+  j)d'  + 

-oo  \ ‘ 

< f)  ' L"  ] 

* «[>  [+G1(*',7'^zi  ^ '^1  ] 


4-  »j< 


i] 


(27) 


Within  the  limits  imposed  by  the  conditions  (12) , Eq. 
(27)  is  an  exact  result.  It  is  not  practical  to  use  it  in 
this  form,  however,  because  the  integrals  must  be  evaluated 
numerically.  In  order  to  reduce  the  computing  time  to  ac- 
ceptable levels,  it  was  found  to  be  necessary  to  perform  at 
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— ...  - - ^ ■»--  >■ 


least  one  additional  analytical  integration.  This  means 
that  approximations  must  be  introduced,  and  consequently  the 
useful  range  of  the  results  is  determined  by  the  accuracy  of 
the  approximations. 

In  order  to  facilitate  the  discussion  of  the  effects  of 
varying  experimental  parameters,  it  is  convenient  to  use  the 
ideal  case  of  collinear,  perfectly-overlapped  beams  as  a 
comparison.  The  generated  field  strength  for  this  ideal 
reference  state,  £$(r)  , can  be  obtained  by  setting  bz  = b, 
and  9C  = ijc  = c = d = 0 in  Eq.  (27)  . Noting  that 

(jx  ip)  = b,  - 1 i ( l - f)i  - G,  fp)  (28) 


and  defining 


Qfp)  - L& l(|°)  G?cp)  + p]  , 


we  find 


y ui. 


r i*»*  i r 

Tf  ^'JeYLp&,Vr)(J<r) 


Eq.  (21)  can  be  used  once  again  to  obtain  an  expression 
which  involves  a single  integration: 

E*  (t)  = iirt0  A e \ tip  r cp)  . < 


where 


vap[-  n ' 

Eq.  (27)  can  then  be  written  in  the  form 

E„<t)  = i.irfto  A'e‘at,'*!)2^pPtoi(»  h[)  , 


where 


5(p)  - i ) Ap*/*’  Qff)  1 e*p  [tfeiCi-cos^)?  H-p) 

+ ii1&jVe>Q  j , »>] 

M1  f » J 

IX  _ j j 

GAfr>|0  ^ &t <i>)Gi  (x'.p)  + i ^2  |°  j 

-a> 

>ey  £ Xp  1*’ 

- if--  ^ 

46iu‘,p)L  £r  ' 


- — 


t -(2-Zp- £c)5.a 9c]  +*ju  | j . 


Since  all  of  the  effects  of  variations  from  the  ideal 
case  are  contained  in  the  function  lip),  Eq.  (35)  serves  as 
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a starting  point  for  systematic  investigations  of  the  ef- 
fects of  changing  the  experimental  parameters.  We  examine 
several  specific  cases  below. 

B.  Special  Cases  of  Nonideality 
1.  Different  beam  sizes 


We  first  consider  the  case  in  which  the  two  laser  beams 
are  collinear,  coplanar,  and  confocal,  but  have  different 
diameters.  The  variable  parameter  is  rt*  = I - bz/b,  , with 
all  other  parameters  having  their  ideal  values.  Under  these 
conditions, 

Gj.  tV.f)  - -l'1  O'f)  - G,  fj»)  -It?,  Aj,  (36) 

and  Eq.  (35)  reduces  to 


1(f)  - 1 (it) 


rili 

rU-e 


c%o«r><xV)-l^i2L±!i- < 


The  integral  in  Eq.  (37)  is  of  the  form  of  Eq.  (21) 
and  thus  can  be  integrated  exactly,  giving 


1(f)-  Q(f)''  [&.'p){G.*<f)  + i*i  f ] 

1 ^ n <**>\ i V l+  u *■  ^ 


t[ 


_t*£  7 

4iif  L6,<p)lbf^)-bl0iij+i^ip]  J . (38) 

2 . Non-conf ocal  beams 

Next  we  treat  the  case  of  collinear,  coplanar  beams 
which  are  not  confocal.  Since  the  crossing  angle  is  as- 
sumed to  be  zero,  the  interfocal  distance  is  df  = Zc  + (\  . 
Here 

(jxL%‘ty)  - b,  -hi  (\~f)  fli(icfd)  r Cr^pM  • (39) 

If  we  make  the  identification  J<:p)  can 

be  written  in  a form  which  is  identical  to  Eq.  (37) , so 
that  the  result  of  the  integration  can  be  written  down  im- 
mediately by  analogy  to  Eq.  (38) : 

“}(f)  - Q<p)  1 [&|tp){Gf(p)+  2.Lcl.f  / + ii-i  p] 

’'i,f['rr7t><3tr,al"li) 

_ dl  (Gty)  nidi 

421p[6,<p){61Vp)t  J * (40) 

3.  Offset  beams 

In  this  section  we  derive  the  function  3(p)  for  col- 
linear, confocal  beams  which  have  the  same  size  and  profile 


. 


. 


but  are  offset  by  a distance  yc  in  a direction  perpendicu- 
lar to  the  direction  of  propagation.  We  have 


= b,-2i  lU-p  - &,*<(»)  , Ml) 

so  that  3(p)  becomes 


(42) 


We  again  see  that  the  integral  in  Eq.  (42)  has  the 
form  of  Eq.  (21) , so  that  the  result  can  be  written  down 
immediately.  Substitution  and  simplification  leads  to  the 
result 

3(p)  - &,if)  Q'f)(fc  + J . (43) 

4 . Crossed  beams 

The  most  complex  single-parameter  variation  is  that 


i 


due  to  the  crossing  angle  $c . Assuming  all  of  the  other 

I 

parameters  to  have  their  ideal  values,  we  have 

GiC<',p)=  b.-UiCI-fOtosSt- 


— P (p)  - Li  9c 

The  expression  for  ^(p)  then  becomes 


(44) 


"3(f)  - i ( ir)  **.» St) z (i-p) 


itj  &Vf) 
t-  — -L a.f)c» 


S -Vt  - -y2. 

Jx'lr.fj-lix'  $c 

f l,1  *' 1 f *l4m  *' 

{ftp)  - li  x'*i*  9t } 


-Cc 


*dX 


_ 'J-l 


41‘p  [6,(p)if(p)-Zi  x'sih  J?6|  +ifc;p] 

f\p)  — Li  X ‘3  3c  J . 


(45) 


The  integral  in  Eq.  (44)  cannot  be  evaluated  exactly 
because  of  its  complexity,  so  it  is  necessary  to  resort  to 
an  approximate  method.  Since  the  crossing  angles  of  experi- 
mental interest  are  typically  small  for  gases,  a natural 
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choice  is  to  make  the  assumptions 


and 


)t‘  f>m  3c 
~ Cl  p) 


JLix'  B| 

Lj « \ 

B,  (f>)  \ r i ti-i  p 


(46a) 


(46b) 


These  conditions  are  obviously  satisfied  for  sufficiently 
small  However,  the  integral  over  'Ji'  covers  the  entire 

axis  (-00,00),  so  that  they  cannot  be  satisfied  uniformly  if 
Sc  ^ 0.  Fortunately,  there  are  other  terms  in  the  integrand 
which  decrease  exponentially  as  /*' | increases,  so  that  the 
conditions  (46)  can  be  assumed  to  hold  in  those  regions 
where  the  integrand  makes  significant  contributions  to  the 
value  of  the  integral.  In  practice,  the  range  of  validity 
of  Eqs.  (46)  can  be  determined  readily  by  a few  simple  nu- 
merical computations. 

When  the  conditions  (46)  are  satisfied,  the  terms  of 
the  form  (|  - oc>'j  ” , where  n = *s  or  1,  in  Eq.  (45)  can  be 
represented  by  truncated  series  expansions.  Numerical  ex- 
perimentation showed  that  for  crossing  angles  within  the 
range  of  experimental  interest  (3C$>  3°),  the  linear  approx- 
imation 

(1  - aC  X')"°  » \ + n<k  %'  (47) 

resulted  in  errors  no  larger  than  about  1%.  This  was  veri- 
fied in  test  cases  by  computing  the  corrections  due  to  the 
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quadratic  terms.  jtroji  JOPY  jrutwjSK®  K> DDG  — - 

Defining 

p lQ(p) 


+ LKtfirK  2 „i 

-4l*p  ‘ty’1*  'Pi  ’ 

(A,  - Zi  5\<n  dc  , 


(48) 


(49) 


and 


r‘f)=  nr,+ fefr  ’ 

we  can  write  Eq.  (45)  in  the  more  compact  form 

i j co  | 

3/p)  = Af  ^i‘p)  Jdy'  |j  (p)  - d X']  [ i /jo)  .(iy| 


(50) 


V £x|0 


4iys7f)^r)-‘,’£'({r'<p)-aW 

_ £t[-y  'Co5^<.  - i^-p)  Sun  ] 


Hp)  - 


ax' 


(51) 


Then,  using  Eq.  (47) , 
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M.^'r  rF",  + IrTriWrY 


ti  \ / ^ ^ t X-  . * *,  V'  I ^ | 

rrr,  K*p( " -r—  -c  ^ ^ *■  ) 


',-•  6,(P)  ti  ,.  4i  f^2  (1 

*£np  i 1 -— *-  - TT-.cc^  C7c  + — 1 — ; — ■*' 

1 1 'f  nf)"  i 


+ \ -l  tjSi«3c.  - — ^ + 1JlilLtl  i.nSi  Ut-Ji 

- *f  r.f) 


£>  ■foi***  / 1 _ 1 lf> 

iix^  rip  l rV> 


__  Zc  ^ lzi\-c)K  h 


r^-  W$t  | *' 


By  expanding  the  last  exponential  term  in  the  integrand 
of  Eq.  (52) , it  is  possible  to  reduce  3ic)  to  the  form 

X 

“ cc 

where  . </,  - ; j,1  * r\.'i 

A.--Afa..r)"'‘[«r,r‘f>i  M'TTjfa^iT,,  (s4) 

- f-j"'f>fV6'p)S,W£  J 


FKOU 


' _L_ 

X,  - l S.rf  7,  ^ pv 


f)  nfj 


\ — - i '-z.  - 

As  — — Dili  5 c £?c 


, - v i r.  -v  x . 2 f; , > ( I - r) 

-l  |.  : -5  >*\  - __- +•  L .J|V1  >L  Cc->  trc 

4r  np> 

_ M b i*\  2c  / i / f-' ' \ 2.i  ft  j.  2 Ci " f- ) . 1u 

ziifr\f)V  r'(fJ  rep)  - 


„ . S.xf  **  ^4it,2« 

^'L-  I — t*-  — -r-. . I.C-5  . «-  “T  — r 

r “j*  r^>‘ 


■j  I rl  ^ V L C-t  ■>  >! 


The  integral  in  Eq.  (53)  can  be  evaluated  by  the  general 


formula 


<c  ti  m 

2_  X IV.  y cu-  ^->1,  x -ut  x z)^* 

x A-  '(T  I 

- Cr.1tV,V|/-.)V!  >» 


,vv  , (Vm-kJ  / ->  . \ yL 

2 u^o 

^ J |\  '. ! (tW  -Kj  ( 

1 1<.  ^ v t n ) 


V,  - 


to  obtain 


Su>  = A , 1 1 - X 


■U|  \ \ / ,*X|  \-  / /(  — 1 / 

. v l5i)  ~ AJ  (fTu)  ( i/i  V /*•  x J • 


5.  Simultaneously  crossed  and  offset  beams 
In  an  actual  experimental  situation,  several  of  the 
parameters  treated  above  may  vary  simultaneously.  To  illus- 
trate the  treatment  of  multivariate  beam  characteristics. 
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we  will  consider  the  case  of  laser  beams  with  a crossing 
angle  9C  and  simultaneous  lateral  offset  ijc  . In  thi 


-s  case 


we  have 


- 1 (it)  f*  4 Q‘f)  2 6-|<)(i  -Cc4»c ) 

+-  i£,'G1%>(4Z1f)"0<fUx24.j:)  ] 

( r ~'/l  I - //. 

* GriCV, jr)  LSiC^)6i(V,p)  + lit  jc  ] 

-oc 


1 .*  L v i a it 

F X -I  til  Sin  cL  — — 

ILF  1 


l X x 


V [S.vpjGiCx'.p;  t * 1Z|0  } 


T.  1 111 

7~~ . j -2(l-p)6«w3t ] (U  H 

(y  1 CX  1 1)/  J J J ) 


(61) 


where  is  again  given  by  Eq.  (44). 

Although  Eq.  (64)  is  rather  complicated,  it  differs 
from  Eq.  (45)  only  in  the  terms  containing  i jc.  in  order  to 
take  advantage  of  this  similarity,  we  write 
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* 


r 


1 i 6,<p)  , 


U f> 


y'  "-i  flx\  $)c  “ 


i% jtx1 

lr 


b 

- ^2  1 * r 
~ t l^c.  (y  1 1 ) | ^ 


(62) 


where  the  last  exponential  term  contains  the  entire  effect 
of  the  offset  >jc*  Making  use  of  Eqs.  (44)  and  (50),  this 
term  can  be  written  in  the  form 


^ <-  v “ * 1'-' 

L EVp)  ( v ) r ‘ p 

r^p)  x ' . 

(63) 
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The  expansion  of  Eq.  (47)  can  now  be  used  again,  reducing 


Eq.  (53)  to 


j 

_ M-1  1 

1 L pip) r 'ipj 

i f3»<p)  r\f>) 

1 

1 x 

1 L_ 

X til 


r r'  c p / I J \ tj 

ft  ‘-0lif>)r<t»r'cf>j  ( Op  + + 1 5^ psp 


^-,1  Jtlllc**'  } . 


i-p)  J *■  ,w  ,J  • (64) 

Therefore,  3tp)  can  be  obtained  directly  from  Eq.  (60)  by 
making  the  replacements 

A, 


r f J 1 P 6a  d 

A|  exf  LU.Cpif^rip)  + 2B,(p)fVp) 


(65) 


and  r'f) 

Lalnpir'ff/ry  + fOMc 


‘fiiffji  ffiVi'i''  1 vf >1 


(66) 


III.  CALCULATION  OF  THE  RELATIVE  POWER  GENERATION 

The  energy  flux  generated  at  a given  point  in  space  is 
proportional  to  the  squared  magnitude  of  the  electric  field 
strength  as  given  by  Eq.  (32)  or  Eq.  (34) . Integrating  over 
the  cross  section  of  the  cell  then  gives  the  total  power 
generation  at  the  frequency  u)^  = 2u),-u>;.  . Alternatively, 
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it  is  possible  to  compute  the  profile  of  the  generated  sig- 
nal over  the  cross  section  of  the  cell  or  to  follow  the 
buildup  of  the  signal  along  the  direction  of  propagation  of 
the  laser  beams. 

The  primary  interest  here  is  in  the  performance  of  the 
system  relative  to  the  ideal  case  (collinear,  confocal  beams 
with  identical  profiles) . This  case  is  somewhat  simplified 
by  the  cancellation  of  the  common  multiplicative  factors  in 
Eqs . (32)  and  (34).  The  major  computational  difficulty  is 
the  multiple  numerical  integration  which  must  be  performed 
in  order  to  obtain  the  relative  power  generation. 

The  numerical  computations  were  performed  on  the  ASC 
system  at  the  Naval  Research  Laboratory.  The  integrals 
were  performed  by  straightforward  application  of  Weddle's 
Rule  (9) . The  only  innovation  was  the  inclusion  of  a vari- 
able step  size  feature  to  increase  the  speed  of  the  integra- 
tion routine  in  regions  where  the  integrand  is  slowly  vary- 
ing. This  modification  was  found  to  be  absolutely  necessary 
in  order  to  obtain  acceptable  accuracy  at  reasonable  cost. 
Typical  computation  times  were  of  the  order  of  30  seconds  of 
CPU  time  for  each  set  of  experimental  parameters.  A more 
detailed  description  of  the  computational  procedure  is  given 
in  the  Appendix. 


, 
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IV.  RESULTS 


A.  Different  Beam  Diameters 

Figure  2 shows  the  results  of  varying  the  parameter 
OCj.  = 1 - bi/bj.  The  negative  values  correspond  to  the  case 
in  which  the  waist  diameter  of  the  lOi  beam  is  larger  than 
that  of  the  cd,  beam,  while  the  positive  values  result  from 
tighter  focussing  of  the  Uiz  beam.  All  of  the  other  experi- 
mental parameters  are  held  constant. 

It  is  quite  clear  from  Figure  2 that  it  is  possible  to 
realize  substantial  relative  power  gains  by  tightening  the 
focus  of  the  uj^  beam  relative  to  the  ib|  beam.  There  are 
experimental  limitations  on  the  extent  to  which  this  can  be 
achieved  in  practice,  but  power  enhancements  of  as  much  as 
50%  (requiring  b;.  ~ 0.65  b,)  seem  quite  reasonable. 

B.  Non-confocal  Beams 

The  variation  in  signal  power  which  results  from  a 
relative  longitudinal  displacement  of  the  foci  of  the  two 
laser  beams  is  illustrated  by  Figure  3.  Positive  values  of 
cl j correspond  to  shifts  of  the  focus  in  the  direction  of 
propagation.  For  the  parameter  values  chosen  here,  the  re- 
sults indicate  that  the  power  generated  by  collinear  beams 
can  be  enhanced  by  about  9%  by  displacing  the  focus  of  the 
USx  beam  by  0.04  cm  in  the  direction  of  propagation. 

C.  Offset  Beams 

Figure  4 shows  the  effect  of  lateral  displacements  of 
two  collinear,  confocal  beams  with  identical  profiles.  The 
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relative  total  power  generation  is  plotted  as  a function  of 
the  offset  parameter  OCj  = ^c/b|#  where  i^c  is  the  actual  off- 
set and  b|  is  the  reduced  beam  radius  defined  by  Eq.  (2). 
Although  the  power  generation  drops  sharply  with  increasing 
0i>,  for  >j  0.01,  the  model  predicts  a slight  increase  in 
power  for  smaller  displacements.  The  maximum  occurs  at 
\j^c  "z.  0.005  b|.  An  important  practical  consequence  of  this 
result  is  the  fact  that  it  is  not  necessary  to  achieve  pre- 
cise collinearity  in  the  alignment  of  an  experimental  CARS 
system. 

D.  Crossed  Beams 

Figure  5 is  a plot  of  the  relative  total  power  genera- 
tion as  a function  of  the  crossing  angle  for  two  non-collin- 
ear  beams  which  cross  at  their  focal  points.  As  expected, 
there  is  a loss  of  power  generation  when  the  beams  are  cros- 
sed. However,  the  magnitude  of  the  decrease  is  not  as  large 

as  might  be  anticipated.  The  model  predicts  that  the  beams 

o 

can  be  crossed  at  angles  as  large  as  ~ 2.5  with  no  more 
than  50%  power  loss.  In  many  cases  this  is  quite  adequate 
to  obtain  a useful  spatial  separation  of  the  CARS  output 
beam. 

E.  Crossed  and  Offset  Beams 

As  an  example  of  the  simultaneous  variation  of  more 
than  one  experimental  parameter,  Figure  6 displays  the  re- 
sults of  variations  of  the  crossing  angle  9«  and  the  offset 
parameter  ^- / b» . it  can  be  seen  that  the  effects  of 
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the  two  parameters  are  essentially  independent  for  small 
variations,  although  there  is  a slight  increase  in  the  opti- 
mum offset  of  the  two  beams  as  the  crossing  angle  increases. 
The  convergence  and  eventual  crossing  of  the  curves  is  of 
little  practical  consequence,  since  the  corresponding  para- 
meter values  are  either  outside  the  region  of  experimental 
interest,  near  the  limits  of  validity  of  the  model,  or  both. 

F.  Longitudinal  Growth  of  the  CARS  Signal 
It  is  easy  to  modify  the  basic  computer  program  to  cal- 
culate the  accumulated  power  over  the  cross  section  of  the 
cell  as  a function  of  the  longitudinal  distance  from  the 
focal  point  in  the  direction  of  propagation  of  the  primary 
beam,  uf,  . The  center  of  the  uJ,  beam  is  assumed  to  propagate 
along  the  centerline  of  the  cell.  This  computation  serves 
to  identify  the  active  volume  of  the  cell  in  which  the  major 
fraction  of  the  signal  power  is  generated. 

Figure  7 shows  the  signal  buildup  starting  from  the 
focal  point  and  proceeding  in  the  direction  of  propagation. 
Results  are  plotted  both  for  collinear  beams  and  for  a cros- 
sing angle  of  0.01  radians.  The  units  of  power  are  arbitrary, 
with  the  reference  level  being  chosen  as  the  cumulative  power 
generated  by  collinear  beams  at  the  focal  point. 

It  is  clear  from  the  figure  that  most  of  the  power  is 


generated  within  a very  narrow  band  (~0.5  cm)  near  the  cen- 
ter of  the  cell.  This  is  a graphic  demonstration  of  the 
increased  spatial  resolution  which  can  be  obtained  by  tight 


focussing  of  the  laser  beams.  The  growth  of  the  signal  is 
roughly  linear  through  the  active  band.  As  might  have  been 
anticipated  from  the  results  shown  in  Figure  3,  the  point  of 
maximum  signal  generation  is  displaced  from  the  focal  point 
in  the  direction  of  propagation. 

Interestingly,  the  model  predicts  a slightly  enhanced 
rate  of  signal  generation  in  the  active  zone  for  the  crossed 
beams  relative  to  the  collinear  case.  The  loss  of  signal 
power  resulting  from  the  beam  crossing  is  thus  attributed  to 
increased  attenuation  as  the  beams  pass  through  the  remain- 
der of  the  sample  cell.  Although  a relative  improvement  in 
spatial  resolution  is  an  expected  result  of  the  beam  cross- 
ing, one  might  expect  to  observe  a noticeable  loss  of  abso- 
lute signal  power  for  an  angle  as  great  as  0.01  radian. 

• G.  Spatial  Profile  of  the  CARS  Beam 

The  basic  computer  program  can  also  be  modified  readily 
to  permit  the  calculation  of  the  generated  radiation  flux  as 
a function  of  the  lateral  position.  This  modified  program 
was  used  to  determine  the  profile  of  the  generated  signal  as 
it  emerges  from  the  exit  face  of  the  sample  cell.  Calcula- 
tions were  performed  for  collinear  beams  as  well  as  for  a 
crossing  angle  of  0.035  radian  (~2°),  and  the  results  pro- 
duced no  surprises.  With  crossed  beams,  the  center  of  the 
CARS  beam  emerged  at  the  point  predicted  by  the  simple  vec- 
tor addition  of  propagation  vectors  and  there  appeared  to 
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be  little  or  no  distortion  of  the  beam  shape.  At  the  end  of 
the  cell  (5  cm  beyond  the  focal  point  of  the  primary  laser 
beams)  the  CARS  beam  was  displaced  from  the  center  of  the  al, 
beam  by  0.08  cm. 


£ ®V.  SUMMARY 

£j  g A computer  model  has  been  developed  to  predict  the 

5 ^ 

ary  relative  signal  power  generated  by  the  CARS  process  under  a 

S3  ' variety  of  experimental  conditions.  It  has  been  demonstra- 

g ted  that  useful  results  can  be  obtained  at  modest  cost  and 

® £ that  the  model  can  be  used  to  aid  the  selection  of  optimum 

2 8 

n a!  experimental  parameters.  In  particular,  the  model  calcula- 
el  tions  indicate  that  precise  alignment  of  the  laser  beams  is 
not  a critical  requirement.  The  predicted  power  losses 
which  result  from  small  lateral  displacements  or  shifts  of 
the  focal  points  are  relatively  modest.  Furthermore,  it 
appears  that  the  laser  beams  can  be  crossed  to  achieve  spa- 
tial filtering  of  the  CARS  signal  without  excessive  signal 
attenuation. 

The  primary  limitation  of  the  model  described  in  this 
paper  is  its  restriction  to  laser  beams  with  Gaussian  pro- 
files. Many  current  CARS  experiments  involve  the  use  of 
beams  with  substantially  different  forms.  For  example,  the 
so-called  "unstable  resonator"  (10)  generates  a laser  output 
with  a cross  section  in  the  shape  of  a doughnut.  The  result 
is  a more  uniform  power  distribution  which  can  produce  a lar- 
ger total  CARS  signal  without  saturating  the  active  medium 
(11) . In  order  to  keep  pace  with  these  experimental  devel- 
opments, the  theoretical  model  must  be  extended  to  laser 
beams  with  arbitrary  cross  sections.  It  is  hoped  that  this 
paper  will  stimulate  further  work  on  the  generation  of 
higher-order  Raman  signals  by  non-Gaussian  laser  beams. 
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APPENDIX 

Details  of  the  Computational  Program 

The  electric  field  strength  E3(f)  of  the  CARS  signal  at 
the  point  r (in  the  coordinate  system  defined  by  Figure  1) 
is  given  by  Eq.  (32)  for  the  ideal  case,  and  by  Eq.  (34)  for 
the  more  general  case.  The  relative  signal  power  generated 
at  the  point  ? is  then  given  by 

“TT  / f \ — - — * ; (A— 1) 


where  P if)  is  defined  by  Eq.  (33)  and  Ilf)  by  Eq.  (35). 

The  total  relative  power  generation  Pre^(z)  at  the 
longitudinal  position  z is  obtained  by  integrating 
over  the  area  of  the  cell: 


The  resulting  three-fold  integral  was  evaluated  numerically 
on  the  ASC  computing  system  at  the  Naval  Research  Laboratory. 
The  FORTRAN  program  was  constructed  in  such  a way  as  to  per- 
mit easy  modification  for  the  purpose  of  computing  trans- 
verse power  profiles,  longitudinal  buildup  of  the  signal, 
etc.  The  program  structure  consists  of  the  main  program 
plus  a set  of  three  nested  subroutines,  as  outlined  below: 
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Main  Program  - initializes  data  and  computes  integral 
over  $ 

Subroutine  RINT  - computes  integral  over  R 
Subroutine  PINT  - computes  integral  over  p 
Subroutine  PGEN  - computes  integrands 

With  this  program  structure,  it  is  very  easy  to  modify 
the  computations  to  deal  with  different  combinations  of  ex- 
perimental parameters.  Since  all  of  the  effects  of  devia- 
tions from  the  ideal  case  are  contained  in  the  factor  sCf) 
in  Eq.  (A-l) , it  is  only  necessary  to  modify  the  subroutine 
PGEN.  One  exception  was  the  computation  of  the  transverse 
profile  of  the  generated  signal,  which  required  minor  alter- 
ations to  the  main  program  and  the  subroutine  RINT. 

Because  of  the  three-fold  numerical  integration,  it 
was  necessary  to  exercise  some  care  in  setting  up  the  inte- 
gration routines  in  order  to  keep  the  computing  time  within 
reasonable  bounds.  The  integrals  were  evaluated  by  Weddle's 
Rule  (9) , which  appeared  to  be  a satisfactory  compromise  be- 
tween maximum  accuracy  and  reasonable  time  requirements.  In 
order  to  keep  the  computational  times  within  acceptable  lim- 
its, it  was  found  to  be  necessary  to  use  a variable  integral 
step  size. 

The  most  effective  strategy  for  minimizing  computation- 
al time  is  to  use  an  automatic  criterion  for  selecting  the 
step  size  in  accordance  with  the  rate  of  change  of  the  inte- 
grand. Although  some  experimentation  was  performed  with 
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such  a procedure,  it  was  found  not  to  be  worth  the  addition- 
al programming  effort  for  the  parameter  combinations  con- 
sidered here.  Instead,  the  integrands  were  plotted  out  for 
a sampling  of  the  variable  parameters,  and  a set  of  inte- 
gral steps  was  selected  to  optimize  the  least  favorable 
case.  The  same  integration  subintervals  and  step  sizes 
were  then  used  for  all  parameter  combinations.  This  pro- 
cedure resulted  in  reductions  of  computing  time  by  as  much 
as  a factor  of  100  in  comparison  with  fixed-step-size  cal- 
culations of  the  same  accuracy.  Since  this  brought  the  to- 
tal computing  requirements  within  acceptable  limits,  no  at- 
tempts were  made  to  achieve  additional  economies.  In  the 
"production"  version  of  the  program,  a single  step  size  was 
used  for  the  integral  over  9,  while  the  integral  over  R was 
divided  into  2 subintervals  and  the  integral  over  the  vari- 
able p was  divided  into  3 subintervals. 
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Fig.  3 — Dependence  of  relative  signal  power  on  the  longitudinal 
displacement  of  focal  points,  df 


Fig-T  — Longitudinal  growth  of  the  CARS  signal  for  collinear  laser 
beams  (o)  and  for  beams  crossed  at  an  angle  of  0.01  radian  (x) 


